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Abstract 

We extend recently proposed variational coupled-cluster method to describe excitation states of 
quantum many-body interacting systems. We discuss, in general terms, both quasiparticle exci- 
tations and quasiparticle-density-wave excitations (collective modes). In application to quantum 
antiferromagnets, we reproduce the well-known spin-wave excitations, i.e. quasiparticle magnons 
of spin ±1. In addition, we obtain new, spin-zero magnon-density-wave excitations which has been 
missing in Anserson's spin-wave theory. Implications of these new collective modes are discussed. 

PACS numbers: 75.10.Jm, 31.15.Dv 
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I. INTRODUCTION 



In the preceding papers (Ref. 1 and 2, hereafter referred to as paper I and II), we proposed 
a general variational theory for ground states of quantum many-body interacting systems. 
Our analysis extends the well-established coupled-cluster method (CCM) {3], \a, Q to a vari- 
ational formalism in which bra and ket states are now hermitian to one another, contrast 
to the traditional CCM where they are not [6|. Ever since the CCM was first proposed, 
attempts have been made to extend it to a standard variational formalism, for examples, 
in the seventies in nuclear physics 7| and later in quantum chemistry [8|. It is perhaps 
fair to say that progress of this variational approach is slow, particularly when comparing 
with a plethora of applications made by the traditional CCM over the last 35 years 
Main difficulties in this variational approach include ad hoc approximation truncations and 
slow convergent numerical results. In I and II, we provided a new systematic scheme to 
overcome these difficulties. In particular, we introduced two sets of important bare distri- 
bution functions and derived self-consistency equations for these functions; calculations of 
physical quantities can all be done in terms of these functions. This strategy is similar to 
that employed by another well-established variational theory, the method of correlated basis 
functions (CBF) 1(|, where density distribution functions are key ingredients. We showed 
that the traditional CCM is a simple linear approximation to one set of bare distribution 
functions. We introduced diagrammatic techniques to calculate those distribution functions 
to high orders for achieving convergent results; resummations of infinite (reducible) diagrams 
can now be done by a practical, self-consistent technique. Furthermore, in our diagrammatic 
approach, a close relation with the CBF method was established and exploited; a possible 
combination of these two methods was also proposed. We demonstrated the efficacy of our 
variational method by appl ying to quantum antiferromagnets. The ground-state properties 
of spin-wave theory (SWT) [ll| was reproduced in a simple approximation. Approximation 
beyond SWT by including higher-order, infinite sets of reducible diagrams produced con- 
vergent, improved numerical results for square and cubic lattices and, interestingly, it also 
cures the divergence by SWT in one dimensional system. 

In this article we extend our variational CCM to describe excitation states. A brief re- 
port of some preliminary results has been published [ijj]. We investigate two different types 
of excitation states using two approaches. In the first approach, we follow the traditional 
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CCM [13I, Q to investigate quasiparticle excitations, but keeping our ket and bra excited 
states hermitian to one another. We then investigate collective modes by adapting Feyn- 
man's excitation theory of phonon-roton spectrum of helium liquid 15] to our method. In 
application to antiferromagnets, we find that quasiparticle excitations correspond to An- 
derson's spin-wave excitations which are often referred to as magnons with spin +1 or —1 



ll| . We find collective modes in these quantum antiferromagnets as longitudinal, spin- zero 



magnon-density-wave excitations which have been missing in Anderson's theory. In our 
approximation, energy spectra of these spin-zero excitations show a large gap for a cubic 
lattice (3D) and are gapless in a square lattice (2D). These spectra are similar to those of 
charge- density- wave excitations (plasmons) in quantum plasmas such as electron gases at 



low temperature 



161 ] . More discussion on these collective modes will be given in the final 



section of this article. 



II. GROUND STATES BY VARIATIONAL COUPLED-CLUSTER METHOD 

We briefly summarize in this section our variational approach for the ground state of a 
many-body interacting system. Details can be found in I and II. We take a spin-s antifer- 
romagnetic Heisenberg model on a bipartite lattice as our model system. The Hamiltonian 
is given by 

l,n 

where index I runs over all lattice sites, n runs over all z nearest-neighbor sites. We use 
Coester representation for both ket and bra ground-states, and write 

|^) = e 5 |$), S = Y J F I C\- <^| = <$|e 5 , S = Y,FjPi, (2) 

1 1 

where model state |$) is given by the classical Neel state, C\ and Cj with nominal index 
I are the so-called configurational creation and destruction operators and are given by, for 
the spin lattice of Eq. (1), 

N/2 - - + + 

£w) = E E h,.,n- il ,f£ jk , (3) 

for the ket state. The bra state operators are given by the corresponding hermitian conjugate 
of Eq. (3), using notation Fj = fi 1 ...j 1 ... for the bra-state coefficients. As before, we have 
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used index % exclusively for the spin-up sublattice of the Neel state and index j for the spin- 
down sublattice. The coefficients {F/, F/} are then determined by the standard variational 
equations as 

The important bare distribution functions, gj = (Cj) and gj = (Cj), can be expressed in 
self-consistency equations as 

9 1 — G(gj, Fj) , g I = G(gj,F J ), (5) 

where G is a function containing up to linear terms in gj (or gj) and finite order terms in 
Fj (or Fj). The Hamiltonian expectation (H) of Eq. (4) can be expressed as, in general, a 
function containing up to linear terms in gj and gj and finite order polynomial in Fj (or in 
Ft), 

{H)=H(g I ,g I ,F I )=H(g I ,g I ,F I ). (6) 

In I and II, as a demonstration, we considered a simple truncation approximation in which 
the correlation operators S and S of Eqs. (2) and (3) retain only the two-spin-flip operators 
as 

The spontaneous magnetization (order parameter) in this two-spin-flip approximation is 
given by the one-body density function as 

(st) = S~p, P = Y.Pij=Y. fijfjij ' ( 8 ) 
3 3 

where we have taken the advantage of translational invariance of the lattice system. For 
j-sublattice, (sp = p — s. Within this approximation the SWT result for the correlation 
coefficient can be derived from Eq. (4) as 

1 



fq — fq — 



V 1 " (^) 2 - l \ » 7, = ;E^ rn , (9) 

J z n 



where f q is the sublattice Fourier transformation of with q restricted to the magnetic 
zone, z is the coordination number of the lattice, and n is the nearest-neighbor index. Fourier 
component of the one-body bare distribution function is derived as, 



Finally, the two-body distribution functions is approximated by, in the same order, 

9ij,i'j' ~ 9ij9i'j' + 9ij'9i'j ■ (H) 

Approximation beond these SWT formulas produced improved results and were given in 
details in II. For simplicity of our first attempt to discuss excitation states, we shall restrict 
ourselves to these approximations of Eqs. (7-11) in the following. 



III. QUASIPARTICLE EXCITATIONS 

As mentioned in Sec. I, inspired by the close relation between our approach and the CBF 
method, we can investigate quasiparticle-density-wave excitations by adapting Feynman's 
excitation theory, as well as usual quasiparticle excitations by similar approach as in the 
traditional CCM. One well-known example of a quantum system exhibiting similar two kind 
of excitations is quantum electron gases [l^ . where quasiparticle excitations are electron 
or hole excitations and collective modes are plasmon excitations representing longitudinal, 
charge-neutral density fluctuations of those quasielectron and holes. In this section we focus 
on quasiparticle excitations and leave discussion of collective modes in the next section. We 
will first discuss these excitations in a general term and then apply to the spin-lattice model 
of Eq. (1) as a demonstration. 



Following Emrich in the traditional CCM 



14|, we express excitation ket-state |\l/ e ) by 



a linear operator X constructed from creation operators acting onto the ground state 



as 



|$ e )=X|^)=Xe s |$), X = ]Tx L Cl, (12) 

L 

and, unlike the traditional CCM, our bra excitation state is the corresponding hermitian 
conjugate, involving only destruction operators as 

(* e | = (* g \X = ($|e*X , X = J £x L C L . (13) 

L 

In Eqs. (12) and (13), the ground-state operators S and S are as given by Eqs. (2), xl 
and its hermitian conjugate xl are excitation coefficients. For quasiparticle creation and 
destruction operators C\ and Cl of Eqs. (12-13), we use index L to mark the following im- 
portant difference to the ground-state counterparts C\ and Cj of Eqs. (2). Due to symmetry 
consideration, some configuration operators are not included in the correlation operator S 



and S of the ground states but they are important in the excited states. In our spin lattice 
example of Eq. (1), the ground-state operators of Eq. (3) always contain even number of 
spin-flip operators (each spin-flip-up operator for the i-sublattice always pairs up with one 
spin-flip-down operator for the j-sublattice) to ensure the total z-component of angular mo- 
mentum s* otal = 0. For the excitation operators, however, the constraints are different. The 
single spin-flip operator for the i-sublattice (or s~j" for the j-sublattice) will be the im- 
portant first term in Eq. (12) to be discussed in the followings; the corresponding excitation 
state |\l/ e ) is in the s* otal = —1 sector (or +1 if s'j' is used). Therefore, these excitations are 
referred to as quasiparticles carrying spin ±1. For our spin lattice models, we expect that 
these quasiparticles are the well-known magnons of spin- wave excitations [11]. 

If the ground state \*f? g ) is exact with energy Eq, the energy difference between excitation 
state of Eqs. (12) and the ground state can be written as, 

(* a \XHX\* g ) _ (* g \X[H,X]\* g ) 

which involves a commutation. In general, \^f g ) is not exact but calculated by approxima- 
tions. For our variational ground states of Eqs. (2), Eq. (14) can be shown to remain valid 
after replacing the exact energy Eq by the variational energy E g which obeys the following 
optimal conditions, 

Bf = <fl ) = (^> = m, (15 , 

9i 9i 

derived from Eqs. (4). 

To prove Eq. (14) after replacing the exact Eq by the variational E g , we first express 
the normalization of excited states of Eqs. (12) and (13) as an expectation value in the 
ground-states of Eqs. (2) as, 

I e = (tt e |* e ) = I g (XX) =I g J2 xl>x l (C v C{) , (16) 

L,L' 

where I g = {^> g \^> g ). We now consider a general linear operator O = 0(C\, Cj) (a polynomial 
of Cj and/or Cj), and write 

0\V g ) = Oe 5 |$) = e s O|$) , (17) 
where the similarity-transformed operator, O = e~ s Oe s = 0{C\,Ci), C\ = C\ and 

Cj = e~ s Oe s = Cj+ [Cj, S] + ^[[Cj, S\, S] + ■ ■ ■ , (18) 



which always terminates for a finite order operator Cj. In each term of such O expansion 
series, by shifting all destruction operators Cj to the right, and using the property C/|$) = 0, 
we conclude that only terms containing constants or only creation operators survive. We 
therefore have a general expression 

0(01^)1%) = 0(ClFj)\* g ), (19) 

where 0{C\, Fj) is a function containing up to linear terms in C\ and finite-order terms in 
Fj. We shall refer Eq. (19) as linear theorem in our variational approach as it is useful for 
general analysis. In fact, the important Eqs. (5) and (6) in Sec. II are two specific application 
of this linear theorem. Therefore we can write, for a special case of Eq. (19), 

C IJ Cl\* a )=Y UL {ti\,F I )\* a ) , (20) 

where Y L ,, L (ClFj) is a function containing up to linear terms in C\ and finite-order terms 
in Fj. Using Eq. (20), Eq. (16) can be written as 

I e = I a Y<XL>Y L , L {g I ,F I )x L . (21) 

L,L' 

Combining with the optimal condition of Eq. (15), it is easy to show 

j ^ g \HXX\^ g )=E g . (22) 

Hence, we obtain similar equation to Eq. (14) for the energy difference, 

6 = U* g \XHX\V g ) -E g = k{X[H,X\) . (23) 

We now apply the above formulas to discuss quasiparticle excitations of spin systems of 
Eq. (1). For simplicity, we consider an approximation in which we retain only single-spin-flip 
operators in X and X of Eqs. (12) and (13), 

X ~ Y, x i s i > * ~ E ^4 , ( 24 ) 

i i 

with coefficients chosen as 

Xi = xM) = » ^ = Hq) = \/f e " q ' ri » ( 25 ) 
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to define a linear momentum q. Such an excited state, |\l/ e ) = X\^f g ), is therefore in the 
sector of sf otaI = — 1 and has a linear momentum q. The normalization integral of Eq. (24) 
is easily calculated as 

j- = (XX) = 2Y^XiXi(sl) + 2s Yl Xi'Xifi'jSlij ~ •''-'•'''/''.; •/'';'//'; > ( 26 ) 

and using Eqs. (8) and (11), we derive, 

£ = 2(a-p)(l + p fl ), (27) 
where p q = f q g q . Using approximations of Eqs. (9-11), we obtain, for isotropic point A — 1, 

hoc-, q^O, (28) 

q 

in all dimensions. 

Calculation of the numerator in Eq. (23) is slightly more complicated. We quote the 
result here as, to the order of (2s) 2 , 

(X[H, X]) w 2s 2 z(l + l Pq + lq9q ) . (29) 

The energy spectrum of Eq. (23) is therefore given by, to the order of (2s), 

e q = ±(X[H,X\) « SZ 1 + 1 P« + ^ . (30) 

Using Eqs. (9-11), we obtain the energy spectrum as 
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« V 1 " (^) 2 ' ( 31 ) 



which agrees exactly with the spin-wave theory 



111 ]. Spectrum of Eq. (31) is gapless in any 



dimension because e q oc q as q — > 0. Similar calculations using spin-flip operators sj and s 
for the j-sublattice in Eq. (24) will produce the same spectrum as Eq. (31) except that the 
corresponding excitation state has spin sf otal = +1. These spin- wave excitations are often 
referred to as magnons. 



IV. QUASIPARTICLE-DENSITY-WAVE EXCITATIONS 

In previous section, by using quasiparticle operators (i.e., spin-flip operators s 1 * 1 ), we have 
reproduced the magnon excitations with spin equal to +1 or — 1. These quasiparticles in 
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general interact with one another, thus producing quasiparticle density fluctuations. Exci- 
tation states due to these fluctuations are usually best discussed in terms of corresponding 
density operator. For our spin models, density operators are clearly given by operators s z as 
it measures the number of spin-flips with respect to the Neel model state and its expectation 
value is the order parameter as given by Eq. (8). For general purpose, we use notation C\ 
for the quasiparticle density operators as opposed to the quasiparticle operators Cf; used 
earlier. The efficiency of using density operators to investigate collective modes of a quan- 



tum interacting system was demonstrated by Feynman 



quantum fluid helium-4 15|, who extended Bijl's theory 17] in a much simpler and clearer 



'or the phonon-roton spectrum of 



fashion. Feynman's excitation formula was also derived by Pines for the plasmon spectrum 



of 3D metals 



16| . The 2D plasmon spectrum first derived by Stern 18j can also be derived 



by using density operator as shown in a PhD thesis [19(. It is interesting to note that both 
the CBF method for the ground state and Feynman's theory for excitation states have been 
successfully applied to fractional quantum Hall effects 2(1 121| . Feynman's excitation theory 
is now often referred to as single- mode approximation 21]. 

Following Feynman, we write the quasiparticle density-wave excitation state as 



\^ e )=X^ g ), X° = ^x L C 



o 

Li 



(32) 



where, as defined earlier, C£ are the quasiparticle density operators. The bra state is given 
by the hermitian conjugate of Eq. (32), (\l/ e | = g \X° . Using the same argument as before 
for the quasiparticle excitation of Eq. (23), we obtain a similar equation for the energy 
difference for our collective modes as, 

n In 



IS 



(X°[H,X }} 



(33) 



where J° = (\l/g|\l/g). We notice that, by definition, density operator C° is a hermitian 
operator, (C°)t = C£. By considering a similar excited state X°\^f g ), it is straightforward 
to derive the following double commutation formula, 



^f- ([X°,[H,X ]]) . 



(34) 



The double commutation in the above equation is the key to the efficiency of Feynman's 
excitation theory. It is often referred to as /-sum rule in other quantum systems such as 



electron gases 



16| 



Before we apply Eq. (34) for collective modes in spin lattices, it is useful to discuss sum 



rules in our spin models as density operators normally obey sum rule equations 
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22] 



The order parameter of Eq. (8) can also be calculated through two-body functions as 

~ (%\%) ' ( ' 

where s z a = l) z sf /N is the staggered spin operator. We introduce total magnon- density 
operator hi as 

2fa = 2s - si + -J2 (36) 

Z n=l 

where as before summation over n is over all z nearest neighbors. Hence, the sum rule for 
the one-body function is simply jrJ2i(ni) — P- The two-body Eq. (35) can now be written 
as, using translational invariant property pi = p, 

2 N/2 

— Y i {n l h i/ )=pp i = p 2 , (37) 

i'=l 



which is the familiar two-body sum rule equation jl0|, l22j . In the approximation of Eq. (7-11) , 
we find that this sum rule is obeyed in both cubic and square lattices in the limit N — > oo. 
In particular, we find that (j| Z)i'(^^i') ~ P 2 ) ^ 1/N in a cubic lattice and oc (In N)/N in a 
square lattice. These asymptotic properties are important in the corresponding excitation 
states as will eb discussed later. However, Eq. (37) is violated in the one-dimensional model, 
showing the deficiency of the two-spin-flip approximation of Eq. (7) for the one-dimensional 
model. We therefore leave further investigation elsewhere and focus on the cubic and square 
lattices in the followings, using approximations of Eqs. (7-11). 

We therefore write our magnon-density-wave excitation state using total magnon density 
operator hi of Eq. (36) as, 

|v|/°>=X 9 |vig, XJ = 5>(?K g>0 (38) 

i 

and its hermitian counterparts for the bra state, (^\ = (^ g \X®. The coefficient, Xi(q) = 
— e iq ri , etc. The condition q > in Eq. (38) ensures the orthogonality between this excited 



state with the ground state. The excitation energy difference is given by Eq. (34) as 



o N(q) 
S°(q) 



^ = ^> ?>0 (39) 
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where N(q) = ([Xj, [H,X°]])/2 and S°(q) = (X° q X° q ) is the structure function. Both N(q) 
and S°(q) can be straightforwardly calculated as, using approximations of Eqs. (7-11), 

sz 

N (<l) = "Y E(7 9 ' + lqlq-<t)9<t, (40) 

and 

= ^(1 + 7> + ^ £[(1 + tJWaW + 27^'^-,'] , (41) 
q' 

where q > 0. The energy spectrum e° of Eq. (39) can then be calculated numerically. We 
notice that Eq. (41) is closely related to the the sum rule Eq. (37) which correspond to q = 
case [with an additional term in Eq. (41) when q — > 0]. Using the approximation of Eqs. (9) 
and (10), it is not difficult to show that N(q) of Eq. (40) has a nonzero, finite value for all 
values of q. Any special feature such as gapless in the spectrum e q therefore comes from the 
structure function of Eq. (41), and hence is determined by the asymptotic behaviors of the 
sum rule Eq. (37) mentioned earlier. 

For a cubic lattice, we plot N(q) and S°(q) for two regions of q in Fig. 1. In Fig. 2 
we plot the corresponding spectra of Eq. (39), together with that of magnon excitations of 
Eq. (31) for comparison. As can be seen from Fig. 2, the spectrum e q has a nonzero gap 
everywhere. The minimum gap is about e q ~ 0.96sz at q = (qo, qo, Qo) with go ~ 0.04?t. 
(This is slightly different to that reported in Ref. 12 where detailed calculations in this region 
had no been were done.) This gap is about the same as the largest magnon energy, e q = sz 
at q = (tt/2, tt/2, tt/2) from Eq. (31). At q = (n/2, ir/2, 7r/2), we have the largest energy 
e° q ps 2.92s,2. This is nearly three magnons' energy at this q. At q = (tt, 0,0), we obtain 
e° q ~ 2.56sz. 

For a square lattice, the structure function S°(q) of Eq. (41) has a logarithmic behavior 
lng as q — > 0. This is not surprising as discussed earlier in the sum rule Eq. (37), where 
occurs the asymptotic behavior (In N)/N as iV — > oo. For small values of q, N(q) approaches 
to a finite value, N(q) ps 0.275sz as q — > 0. The corresponding energy spectrum of Eq. (39) 
is therefore gapless as q — > 0. Similar to the cubic lattice, we plot N(q) and S°(q) of a square 
lattice in Fig. 3 and the corresponding spectra of Eq. (39) and Eq. (31) in Fig. 4. As can 
be seen from Fig. 4, magnon- density- wave energy is always larger than the corresponding 
magnon energy. At small values of q (q < 0.057r), we find a good approximation by numerical 
calculations for the structure function, S°(q) ps 0.31 — 0.16 lng with q x = q y . Similar behavior 
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holds near q = (ir,ir). The energy spectrum of Eq. (39) in these region can therefore be 
approximated by 

n 0.275sz . . 

e° « — , q -> 42 

q 0.31-0.161ng' H V ; 

for a square lattice with q x = q y . We notice the slight difference for the coefficients of Eq. (42) 

to that of Eq. (19) of Ref. 12 where we focused in the region with q y = 0. Although our 

calcuations clearly show this spectrum of a square lattice is gapless at q = and q = (ir, it), 

it is nevertheless very "hard" when comparing with the magnon's soft mode e q oc q at small 

q. For example, we consider a system with lattice size of N = 10 10 , the smallest value 

for q is about q ~ 10~ 10 7r and we have energy e° q ~ 0.07sz. Comparing this value with 

the corresponding magnon energy e q ~ 10~ 10 sz, we conclude that the energy spectrum of 

Eq. (39) is "nearly gapped" in a square lattice. We also notice that the largest energy in 

a square lattice e q ~ 2.79sz at q = (ir, 0), not at q = (tt/2, vt/2) as the case in a cubic 

lattice. At q = (n/2, 7r/2), we obtain e° ~ 2.62sz for the square lattice. We will discuss 

physical implications of these excitations in the next section. 



V. DISCUSSION 



We have obtained in this article two main results. Firstly, we have succeeded in extending 
our recently proposed variational approach to describe, in general terms, excitation states 
of a quantum many-body system. Secondly, we have applied our technique to quantum 
antiferromagnets thus reproducing the well-known magnon excitations and, in addition, 
we have obtained a new, spin-zero longitudinal collective modes which have been missing 
in spin- wave theory of Anderson [11]. In the followings, we shall discuss further physical 
implications of these new excitations and we conclude this article with a summary in the 
end. 

It is interesting to notice similar behaviors between collective modes of quantum antifer- 
romagnets and plasmon excitations of electron gases as both spectra show a large energy gap 
in 3D and are gapless in 2D. In fact, further similarity between these two quantum systems 
can be made. It is generally accepted that, for many purposes, a quantum antiferromag- 
net at zero temperature can be considered as a gas of weakly interacting, equal numbers 
of spin ±1 magnons (the transverse spin-flip wave excitations with respect to the classical 
Neel state); also present in the system are the spin-zero, longitudinal fluctuations consist- 
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ing of multi-magnon continuum 23j, [24], |25|, |26|. This is similar to quantum electron gases 



which can also be considered as a gas of weakly interacting, equal numbers of quasielectrons 
and holes (the transverse excitations near the Fermi surfaces) and the charge-neutral, Ion- 

n 

gitudinal fluctuations producing quasielectron-hole continuum [16j. Plasmon excitation of 
electron gases have been well observed as sharp peaks over the electron-hole continuum [if]] . 
However, plasmon-like collective modes of quantum antiferromagnets as discussed in this 
article have so far eluded from observation to our best knowledge. We can only draw some 
support by considering a finite-size Heisenberg model of Eq. (1). As the ground state of a fi- 
nite antiferromagnetic Heisenberg lattice is spin-singlet, we expects lowing- lying excitations 
are triplet with ^-component of spin equal to 0, ±1. As lattice size increases from finite to 
infinite, for the cubic and square lattices, the spontaneous-symmetry-breaking occurs and 
the ground state is no longer spin singlet but has a long-ranged antiferromagnetic order. 
We expect that the triplet excitation splits into different branches. The magnon spectrum 
of Eq. (31) with spin ±1 and the spectrum of Eq. (39) for spin- zero magnon- density waves 
are our approximation for these different branches of excitations. We also notice that re- 
cently modified spin-wave theories were applied to finite systems with results in reasonable 



agreements with exact finite-size calculations 



23, 



28. 



29j. As pointed out in Ref. 28, how- 



ever, a major deficiency in this theory is the missing spin-zero excitations as the low-lying 
excitations for a finite lattice Heisenberg model are always triplet as mentioned earlier. We 
believe our magnon-density-wave excitation as discussed here corresponds to the missing 
branch; the energy gap in the cubic lattice and the nearly gapped spectrum in the square 
lattice of Eq. (39) reflect the nature of long-ranged Neel order in the ground states of infinite 
systems. Improvement for spectra of Eq. (39) can be done in similar fashion as was done 
for the ground state detailed in paper II, particularly for the square lattice. We will have 
more motivation to do so if we have experimental evidence of these collective modes. 

In any case, this artcle concludes our general presentation of a new formalism of vari- 
ational coupled-cluster method for a quantum many-body system. Beginning in paper I, 
we introduced and discussed bare distribution functions, key ingredient of this formalism. 
In paper II, we developed diagrammatic techniques for practical, high-order calculations of 
these functions. Application to quantum antiferromagnets has demonstrated the efficacy 
of this technique. Present artcle extends this formalism to excitation states. As discussed 
earlier, application to quantum antiferromagnets have produced new modes which have been 
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missing in all spin- wave theories and are yet to be confirmed by experiment. Our next main 
focus is to combine our present variational approach with the CBF method as first discussed 
in paper II. Hence we write our new ground state as 

|M/ u )=e s > 9 ) = e 5 V|$), (43) 

where S is as given by Eqs. (2) and S° is the generalized Jastrow correlation operator 
involving quasiparticle density operators as 

S° = Ef^h ( 44 ) 

ij 

with /P- as new variational functions. Using the 2-spin-flip approximation of Eq. (7) for 
S, the new wave function of Eq. (46) can be understood as including both quasipartcle 
fluctuations described by operator exp(S) and quasiparticle-density fluctuations described 
by operator exp(5 <0 ). The results of collective modes obtained in Sec. IV certainly make this 
combination of Eq. (43) much more appealing and imperative. 
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Fig. 1 N(q) and S°(q) of Eqs. (40) and (41) for a cubic lattice. Shown are the values for two 

regions q = (0, 0, 0) to (tt, 0, 0) and to (tt, tt, tt). 
Fig. 2 Excitation energy spectra in unit of sz in a cubic lattice. The higher branch is for the 
plasmon-like excitation of Eq. (39) and the lower one is for the magnon excitation of Eq. (31). 
Fig. 3 Similar to Fig. 1 but for a square lattice. The divergence of S°(q) at q = (0,0) and (tt,tt) 

is given in the text. 

Fig. 4 Similar to Fig. 2 but for a square lattice. The behavior near q = (0, 0) and (it, tt) for 

magnon-density- waves is given by Eq. (42). 
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